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Recently, a framework for controller design of sampled-data nonlinear systems via 
their approximate discrete-time models has been proposed in the literature. In this paper 
we develop novel tools that can be used within this framework and that are very useful 
for tracking problems. In particular, results for stability analysis of parameterized time- 
varying discrete-time cascaded systems are given. This class of models arises naturally 
when one uses an approximate discrete-time model to design a stabilizing or tracking 

■ controller for a sampled-data plant. While some of our results parallel their continuous- 
time counterparts, the stability properties that are considered, the conditions that are 

■ imposed and the the proof techniques that are used are tailored for approximate discrete- 
C*~) ■ time systems and are technically different from those in the continuous-time context. We 

illustrate the utility of our results in the case study of the tracking control of a mobile 
i -^h ' robot. This application is fairly illustrative of the technical differences and obstacles 

d • encountered in the analysis of discrete-time parameterized systems. 

s 

> 

1 Introduction 

& ■ 

The prevalence of digitally controlled systems and the fact that the nonlinearities in the plant 
model can often not be neglected, strongly motivate the area of nonlinear sampled-data systems. 
A typical nonlinear sampled-data system consists of a nonlinear continuous-time plant and 
a nonlinear discrete-time controller that are interconnected via the analog-to-digital (A-D) 
and digital-to-analog (D-A) converters. Despite the importance of this class of systems, few 
systematic tools for nonlinear sampled-data controller design are available in the literature. 
Instead, in control practice one typically follows the commonly accepted hypothesis that if a 
continuous-time controller designed for a continuous-time plant, is implemented at a sufficiently 
fast sampling rate, the sampled-data system should "behave well". While this intuition is 
correct in general, the required sampling may be too fast to be implemented in practice because 
of the available hardware limitations. 
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The stumbling block that naturally arises in the formal analysis and controller design for 
nonlinear sampled-data systems is the fact that the model of the system is rather complex 
(hybrid, nonlinear, periodically time- varying) and very hard to deal with directly. Consequently, 
there are several different methods that can be used when designing the controller. 

One method, which is sometimes referred to as the emulation method, consists of designing 
a continuous-time controller based on the continuous-time plant model and then discretizing 
the controller for digital implementation. We stress that the emulation method ignores the 
sampling during the controller design step. Results on emulation for nonlinear sampled-data 
systems can be found in ^Hj and references defined therein. Another method, which for ease 
of reference we will call here, exact (respectively approximate) discrete-time design (DTD), 
consists of obtaining the exact (respectively, approximate) discrete-time model of the plant 
and then designing the controller based on the discrete-time plant model. In this method, the 
inter-sample behaviour is ignored during the controller design. 

Since the emulation method does not take the sampling into account in the controller design 
step, it is reasonable to expect that DTD method may produce better results than the emulation. 
Unfortunately, if we want to implement the second method, even if the nonlinear continuous- 
time plant model is known, in general we will be unable to obtain the exact discrete-time 
model of the plant since this involves solving a nonlinear differential equation analytically over 
one sampling interval. Instead, we need to exploit an approximate discrete-time model that 
is obtained using some numerical integration scheme, such as Runge-Kutta (see e.g. [HI USD- 
Early results that use approximate discrete-time plant models for controller design can be found 
in jU |H] • See also the more recent works of [23 123 E3 • 

The main pitfall of the approximate DTD method is that if one is not careful with the 
choice of the approximate model and the design of the controller, it is possible that a controller 
asymptotically stabilizes the approximate plant model but not the exact model. It is noteworthy 
that this fact concerns even linear systems as we illustrate next. 

Example 1 Iffffl Consider the system 

%1 — X 2 / -i \ 

x 2 = u 

whose exact and Euler models are respectively 

Euler: I ^ + ]\ = ^S+^ff (2) 

\ x 2 (k + 1) = x 2 (k) + Tu(k) v ; 

J zi(fc + l) = x 1 (k)+Tx 2 (k) + 0.5T 2 u(k) ,„s 
* XaCl - \ x 2 (k + l) = x 2 {k)+Tu{k). {6> 

If we define the controller Ut{x) = — (1/T)[x\ + 2x 2 ] for the Euler model, we get that one 
eigen-value of the closed loop system with the exact-model is on the unitary circle for all values 
ofT > 0. Hence, the system cannot be asymptotically stable for any T > 0. On the other hand, 
the eigenvalues of the approximate model are ±a/ (1 — T) and one can actually show that there 
exists b>0 such that \x(k)\ < b\x(0) \ e~°- 5kT for all x(0) G M 2 and all T E (0,0.5). 

Motivated by this fact, a framework for nonlinear sampled-data controller design via ap- 
proximate discrete-time models has been proposed in [HBl EH ESI EH ED] • These results are 
very similar in spirit to results from the numerical analysis literature (see e.g. jlSl HI) that 
applies to continuous-time control systems. In [33] checkable conditions on the continuous-time 
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plant model, the controller and the approximate discrete-time model are presented which guar- 
antee that if the controller stabilizes the approximate model, it would also stabilize the exact 
discrete-time model. Furthermore, in it was shown that stability of the exact discrete-time 
model under mild conditions guarantees also stability of the sampled-data system. Hence, the 
results of [SHI El] provide a framework for controller design of sampled-data nonlinear systems 
via their approximate discrete-time models. 

The above mentioned results are primarily targeted at establishing conditions for stability 
that is, they are prescriptive and non- constructive. In other words, they provide a framework 
for controller design but without explicit recipes for controller design. A range of different 
constructive methods for controller design within the above given framework has been reported 
in the literature: backstepping via the Euler model of strict feedback systems j2H]j Lyapunov 
methods based on changes of supply rates for input-to-state stable systems jS2 UH UH|; and 
optimization based stabilization [Hj. 

Other constructive design methods for non-parameterized systems can be found, for instance 
in [23J ESI EJ E] , in the survey j2B] and in the references listed therein. 

Contributions of this paper 

The results that we present in this paper contribute to what we may call cascades- based control. 
Roughly speaking, this approach aims at designing controllers in cases when the closed loop 
system has a cascaded structure. Moreover, the closed loop dynamics shall verify certain 
structural conditions imposed either on the functions that define the closed loop dynamics or 
indirectly, in terms of properties of Lyapunov-like functions. While there is a wide number 
of such results establishing different forms of asymptotic stability for continuous-time systems 
(both, autonomous and non-autonomous), there are only a few results for discrete-time systems. 
In particular, we cite ^2] where sufficient conditions for stability of cascades that use the input- 
to-state stability (ISS) property are presented. 

Yet, the study of stability of cascaded systems in the discrete-time context has a double 
motivation: Firstly, it obviously inherits the motivations from the continuous-time case: there is 
a wide range of applications in control design, of stability results for cascades. See 1201 1221 EH] 
and references therein for a large number of results and applications in this area; from a control 
theory viewpoint, the motivation for this research originated probably in geometric nonlinear 
control where it was shown that many systems can be transformed into a cascade via a local 
change of coordinates (see, for example, [TTH Lemma 1.6.1]). Secondly, in view of the previous 
discussions it is highly desirable to establish conditions tailored specifically for approximate 
discrete-time systems, under which one can rely on the DTD control design method. 

This paper is aimed in that direction. We present results on uniform asymptotic practical 
stability of parameterized (in the sampling period) discrete-time cascaded systems. Thereby, 
contributing to the framework established in (221 HH3 El] and other above-cited references for 
sampled-data systems. Although we have been inspired by continuous-time results on stability 
of cascades (in particular by those in [21 EH]), the results presented here are not a simple 
translation of their counterparts in continuous-time. Indeed, the properties we consider, the 
conditions we impose and the proofs we establish here are notably different from continuous- 
time ones and they are tailored specifically for discrete-time parameterized systems that arise 
within the above mentioned framework. 

Furthermore, advantages of DTD control design are illustrated with the sampled-data track- 
ing control problem of the unicycle benchmark problem where we show an improvement in per- 
formance with respect to continuous-time based designs. In particular, we will illustrate how 
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to apply our main results on stability of cascades. While the benchmark problem is the same 
as in continuous time, this case-study is fairly representative of the large differences between 
discrete-time and continuous-time based designs as well as the tools involved in the study of 
the proper stability properties for parameterized systems. Indeed, the proofs of the results 
for the unicycle are original and also very different from their continuous counter-parts. Thus, 
through this case-study we illustrate how using our results we can obtain new control algorithms 
that can be regarded as a continuous-time controller redesign that outperforms the emulated 
controllers. 

The rest of the paper is organized as follows. In Section |2] we provide mathematical prelim- 
inaries with some definitions and basic results for parameterized systems. Our main results are 
presented in Section El in the form of two different theorems that are "Lyapunov-based" and 
"trajectory based" respectively. In Section 0] we present the unicycle case study and the proofs 
of the main theorems are presented in Section E] Some concluding remarks are presented in 
Section Further technical proofs and other auxiliary results are included in the appendix. 

Notation. A function a : M>o — > M>o is said to be of class K. (a G fC), if it is continuous, 
strictly increasing and zero at zero; a G /Coo if, in addition, it is unbounded. A function 
p : M> x M> -> M> is of class KC if for all t > 0, p(-,t) G /C, for all s > 0, (3(s, •) is 
decreasing to zero. A function 7 : IR> — > M>o is said to be of class M if 7(-) is continuous and 
nondecreasing. We denote by |-| the Euclidean norm of vectors. We denote by M and N the 
sets of the real and natural numbers respectively. For an arbitrary r G M we use the notation 
|_rj := max z. Given strictly positive real numbers L,T we use the following notation: 

zGZ,z<r 

(4) 

2 Parameterized discrete-time systems 

In this section we present a result that provides a framework for controller design for sampled- 
data systems via their approximate discrete-time models. This result generalizes Theorem 
1] in that it is applicable to time- varying plants and hence more suitable for tracking problems. 
The proof of this result follows the same steps as that of Theorem 1 in the above paper and 
hence it is omitted. Our main results that are presented in Section 3 facilitate sampled-data 
controller design within the framework of Theorem 1 in this section. 

We also introduce here the precise definitions of stability that we pursue and other prelim- 
inary important results. 

2.1 Controller design framework 

Consider the class of systems 

x(t) = f(t,x(t),u(t)) 

y(t) = h(x(t)) (5) 

where x G M. nx and u G M m are respectively the state and control input. We assume that 
for any given x and u(-) the differential equation in (jSJ) has a unique solution defined on its 
maximal interval of existence [0,t max ). This may be guaranteed, for instance, by requiring / 
in (jSJ) to be locally Lipschitz, uniformly in t. The control is taken to be a piecewise constant 
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signal u(t) = u(kT) =: u{k), Wt G [kT, {k + 1)T), k G N, where T > is the sampling period. 
Assume that some combination (output) or all of the states (x(k) := x(kT)) are available at 
sampling instant kT, k G N. The exact discrete-time model for the plant (JHJ, which describes 
the plant behavior at sampling instants kT, is obtained by integrating the initial value problem 

x{t) = f{t,x{t),u{k)) , (6) 

with given u(k) and xq = x(k), over the sampling interval [kT, (k + 1)T]. 

If we denote by x(t) the solution of the initial value problem (jHJ) at time t with given 
Xo = x(k) and u(k), then the exact discrete-time model of (jSJ can be written as: 

r(k+l)T 

x(k + 1) = x(k) + / f(r, x(r),u(k))dr =: F^(x(k),u(k)) . (7) 

We emphasize that is not known in most cases. Indeed, in order to compute F£ we have 
to solve the initial value problem © analytically and this is usually impossible since / in (JHJ) 
is nonlinear. Hence, we will use an approximate discrete-time model of the plant to design a 
discrete-time controller for the original plant (JHJ) • 

Different approximate discrete-time models can be obtained using different methods, such 
as a classical Runge-Kutta numerical integration scheme for the initial value problem (JBJ) -see 
e.g. [H 123 12H HS1 E] • The approximate discrete-time model can be written as 

x(k + 1) = F«(k,x{k),u{k)) . (8) 

For instance, if / is locally Lipschitz in t and x, the Euler approximate model can be defined 
as x(k + 1) = x{k) + T f{kT,x{k),u{k)) and it can be shown to be an 0(T 2 ) approximation 
of the exact discrete-time model. On the other hand, if / is measurable in t, then a modified 
"Euler" model that is 0(T 2 ) approximation of the exact model is given by x(k + 1) = x{k) + 

f£ +1)T f(r,x(k),u(k))dr. 

In our work, the sampling period T is assumed to be a design parameter which can be 
arbitrarily assigned. Since we are dealing with a family of approximate discrete-time models 
F£, parameterized by T, in order to achieve a certain objective we need in general to obtain 
a family of controllers, parameterized by T. Thus, we consider a family of dynamic feedback 
controllers 

z(k+l) = G T {k,x{k),z{k)) , g , 
u(k) = ur(k, x(k), z(k)) , 

where z G M. nz . We denote the right hand sides of the closed-loop systems (JHJ), © an d J7J), © 
respectively as 

~ a( _ f F$(k,x,u T (k,x,z)) \ ~ ef . / F*(k,x,u T (k,x,z)) \ 
^ Tlfc ' X) - V G T (k, x,z) )> ^ Tlfc ' X) - \ G T (k, x, z) ) 

where x = (x T z T ) T . 

We emphasize again that if the controller @ stabilizes the approximate model (jHJ) for all 
small T, this does not guarantee that the same controller would approximately stabilize the 
exact model (JU) for all small T. In addition to Example several inherently different examples 
that illustrate this phenomenon can be found in \§\ ITH]. The following result that was 

proved in [HS] for time-invariant systems, gives sufficient conditions for stabilization of (J2J) via 
controllers Q that are designed using (JHJ). The result states that one can conclude semiglobal 
uniform asymptotic stability (SP-UAS) of the exact discrete-time closed loop system from the 
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same property for the approximate discrete-time closed loop and an appropriate consistency 
property between the exact and approximate models (see Definition [TJ. We emphasize that 
under very weak conditions this guarantees that the inter-sample behaviour of the sample-data 
closed- loop system would also be bounded (see [M]). 

The following result generalizes |3*3*l Theorem 1] to time-varying systems. The proof is 
omitted since it follows from that of Theorem 1 in the above paper with straightforward changes. 



Theorem 1 Suppose that there exists (3 G ICC such that for any strictly positive numbers 
[L,r], A, 5) there exist a : M>o x M> — > K>o U oo and T* > such that for all k Q > 0, 
\x(k )\ < A andT G (0,T*) we have' 



SP-UAS of approximate: 

The solutions of (OJ) satisfy: 



\4> a T (k,k ,x(k j)\ < P{\x{k )\,T{k-k )) + 6, 



k > k > 



(10) 



multiple-step consistency between Fj, and F^: 

max.{xi, x 2 } < A ; \x\ — X2I < c and k > imply 



where a satisfies 



F^{k, Xl ) -F£{k,x 2 ) <a{c,T) 



a k (0,T)< V , \/ke[0,£ L<T ] 



(12) 



Then, for any strictly positive real numbers (A, 5) there exists T > such that for all k > 0, 
|x(fc )| < A and T G (0, T) the solutions of (0) satisfy: 

• SP-UAS of exact: 



(k,k ,x(k ))\ < (3{\x{k )\,T{k - ko)) + 5, 



Wk>k o >0. 



(13) 



□ 



Checking SP-UAS for the approximate model (see the equation (|l()jl ) is very hard in general 
and in this paper we provide new results that facilitate checking this property in the case when 
the approximate model has a cascaded structure. 

Remark 1 The consistency condition in Theorem Q is checkable although F^ is not known 
in general. For instance, it was shown that this condition holds if for any A > there exist 
K, T* > and p G K, such that max{|5i| , \xi\} < A, k > and T G (0, T*) imply 



- ^(k,x 2 ) < (1 + KT) \x 1 -x 2 \+Tp{T) 



This condition is readily checkable and it is commonly used in numerical analysis literature \45^ - 
For instance, it holds if the following two inequalities hold in a semiglobal practical sense: 



F^{k,xx) -F£{k,x 2 ) < {l + KT)\xi-x 2 \ 
F^(k,x) - F%{k,x) < Tp{T) . 



(14) 
(15) 
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The condition fli^| ) is a particular local Lipschitz property of F^, which is usually not hard 
to check (see Remark^. The condition 4i<5j) is ensured by choosing an appropriate consistent 
approximate model F£ for controller design. A range of such approximate models can be found in 
standard books on numerical analysis \4$j (if f is independent oft),inJ^/(iff is differentiable 
in t) and in JEj (if f is measurable in t in the Lebesgue sense). 

Theorem and the discussion in the Introduction motivate the following stability property. 
Note that below we use the "max" instead of the "+" characterization of the practical stability 
that was used in Theorem ^ since the "max" formulation is easier to use in our proofs, but the 
two characterizations are qualitatively the same. 

Definition 1 (SP-UAS) The family of the parameterized time-varying systems 

y(k + l) = F T (k,y(k)) (16) 

is semiglobally practically uniformly asymptotically stable SP-UAS (resp. uniformly globally 
asymptotically stable UGAS) if there exists [3 G KL such that for any pair of strictly positive 
real numbers (A, v) there exists T* > (resp. there exists T* > 0) such that for all k a > 0, 
y(k a ) = y Q with \y \ < A, T G (0,T*) (y(k Q ) = y a with y Q G W 1 , T G (0,T*)J the following 
holds: 

\(j)^(k,k ,y )\ < max{f3(\y \ , (k - k Q )T),u} (resp. \4> y T (k, k ,y )\ < f3(\y \ , (k - k Q )T) ) 

(17) 

for all k > k a . □ 

Notice that the convergence of solutions in (JT7j) is not uniform in the sampling period since the 
solutions are allowed to converge slower as T increases. However, the bound (|17|) has uniform 
overshoots. The motivation to consider this precise type of asymptotic stability is twofold. 
First, this is the same property that the approximate model should satisfy in the statement of 
Theorem ^ Second, under the conditions of this theorem, the exact model would also satisfy 
the stability property of Definition [T] and by further applying results of jMj we can conclude 
that the sampled-data systems is SP-UAS. 



2.2 Problem setting and other definitions 

This paper focusses on the problem of establishing sufficient and necessary conditions for SP- 
UAS of cascaded parameterized systems of the general form 

x(k + l) = f T (k,x(k),z(k)) (18) 
z(k + l) = g T (k,z(k)) , (19) 

where x G M. nx , z G M nz and T is the sampling period. For the cascade (|18|). ()19j) we use 
the notation £ := [x T z T ] T to denote the state of the overall system. These results are used 
within the framework provided by Theorem to check appropriate stability properties of the 
approximate discrete-time model. 

To state our main results we also consider the system 

x(k + l) = f T (k,x(k),0) , (20) 
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and often, we will regard z in the system (|18|) as an exogenous input that is not necessarily 
generated by the subsystem (JTHJ). In that case, we refer to the subsystem ffTHjl as the system 
with input z. The solution of the system (j!8j) with input z at time k that starts at initial 
time instant k from the initial state x(k ) = x and under the action of the input sequence 
u? ko fe N := {z(k ), . . . z{k — 1)} is denoted as <^r(/c, k Q , x , u)^ ka k ^). We also use u z := uj? ko s. Note 
that the solution of the system (J20|) is the same as the solution for system (JTBJ) with input z 
when z(j) = 0, Vj G [k a ,k] and hence for solutions of (J2DJ) we use the notation 4>j>(k, k , x a , 0). 
Similarly we use the notation 4>%{k, k , £ ), 0t(^5 £°) an d 0t(^; ^o, z a ) to denote solutions of 
the overall system JUJ), an d its x and z components respectively 

Roughly speaking, our main results establish SP-UAS of the cascade under the condition 
that both subsystems, fl~9|) and (|20|) are SP-UAS and that the solutions of the cascade be 
uniformly bounded in the following specific manner. 

Definition 2 The system A16}) is uniformly semiglobally bounded USB (uniformly globally bounded 
UGB), if there exist k G /Coo o,nd c, such that for any A > there exists T* > (there exists 
T* > 0) such that k > 0, y(k a ) = y Q with \y a \ < A and T G (0,T*) (y a G W l and T G (0,T*)) 
implies 

\<P v T (k,k 0} y )\<K(\y \) + c } (21) 
for all k > k a . □ 

We also use the following property which is very related to SP-UAS. (cf. remark |2J). 

Definition 3 The family of the parameterized time-varying systems Mb\) is Lyapunov SP- UAS 
if there exist a,\, ot2 G /Coo, «3 G /C ; L G M and for each pair (A, u) > there exists T* > 
and for each T G (0, T*) a continuous function '■ K> x IR n — > IR> such that for all \y\ < A, 
all k > and all T G (0, T*) we have that 

axQyl) <V T (k,y)<a 2 (\y\) (22) 

V T (k + 1, F T (k, y)) - V T (k, y) < -T(a 3 (\y\) + v) (23) 

and 

| V T (k, r) - V T (k, s)\< L (max{|r| , \s\}) \r - s\ . (24) 

for a//max{|r| , \s\} < A, T G (0, T*) and k > 0. T/ie system is Lyapunov UGAS if there exists 
T* > snc/i t/icrf t/ie a6ofe conditions hold for all x, r, s G M n ; > and wzt/i i/ = 0. □ 

The Lipschitz condition (|24|) holds in particular if Vr(/c, •) is continuously different iable and the 
derivative of Vr(k, •) is bounded on compact sets, uniformly in small T and k > 0. 

Remark 2 We note that the properties in Definitions^ and\^ are very related. In particular, 
it was shown in W% that 1%%)) . are equivalent to the bound (fT7| ). However, the converse 
Lyapunov theorem in Iffffl does not produce a Lyapunov function satisfying the condition \2$ 
and we believe that constructing such converse Lyapunov functions is an open problem in the 
literature. In particular, it would be important to provide conditions under which one can 
construct Lyapunov functions from Definition^ for families of systems Mb\) with a discontinuous 
right hand side. Such converse theorems for non-parameterized systems can be found in mH 
but we are not aware of similar results for parameterized systems. 

The following property is similar to continuity of solutions of differential equations satisfying 
the local Lipschitz condition if we think of t := kT as "continuous time". This will be crucial 
in establishing a trajectory-based proof of our main result. 
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Definition 4 The solutions of the system hltyl with input z are uniformly semiglobally contin- 
uous USC (uniformly globally continuous UGC) if for any A there exists T* > (there exists 
T* > 0) such that for any r\ G (0, A), e > and L > there exists \i > such that for all 
T G (0, T*), all z(-) with \lv z \ < \i, k D > and all x(k D ) = x with \x \ < r\ we have that: 

|0r(fc, k ,x ,uf kotk) ) - </%(k,k o ,x o ,0)\ < e , VA; G [k a , k a +#l,t\ ■ (25) 

□ 

The last technical tool that we introduce is a set of checkable conditions for USC inspired by 
the literature on numerical methods. The proof can be found in the Appendix. 

Proposition 1 Suppose that the system is SP-UAS (UGAS) and for any pair of strictly 
positive (Ai, A 2 ) there exist strictly positive K,T* such that: For all k > 0, max{|xi| , |x 2 |} < 
A 1; maxjl^xl , \z 2 \} < A 2 and T G (0,T*) we have: 

\f T {k,xx,Zx) - f T (k,x 2 ,z 1 )\ < (l + KT)\xx-x 2 \ (26) 
\f T (k,xx,Zx) - fT(k,xx,z 2 )\ < KT\z x -z 2 \. (27) 

Then, the system \18) with input z is USC (UGC). □ 

Note the slight difference between conditions in Proposition ^ and conditions presented in 
Remark ^ 

3 SP-UAS of time-varying cascades 

We present next two main results that deduce stability of the cascade (JTBJ), (JTHJ) from stability 
of two lower dimensional auxiliary systems (J2*U|) and ()19j) and the boundedness of trajectories 
of the overall system ()18j). f)19|) . These results are fundamental in that they are necessary and 
sufficient for stability of the cascade under appropriate conditions. Hence, the first condition of 
Theorem ^ can be checked via by three separate conditions that are usually easier to verify. Our 
two main results are similar in spirit but they are derived under slightly different conditions. 

Theorem El states that under the USC assumption, we have that USB (UGB) of the system 
(HEJ), (HH and SP-UAS (UGAS) of the auxiliary systems flUJ) and (O is equivalent to SP-UAS 
(UGAS) of the overall cascade (J2"U|) and (fTT?|) . The proof of this theorem is inspired by the 
trajectory based proof in [J2j Theorem 1] that does not appeal to converse Lyapunov theorems 
(see Section EJ). Theorem El uses Assumption d (see below) instead of the USC assumption and 
it shows that USB (UGB) of the system ijIB} . ipj) . Lyapunov SP-UAS (UGAS) of the auxiliary 
system (f^OJ) and SP-UAS (UGAS) of the auxiliary system (HHJ) implies SP-UAS (UGAS) of the 
overall cascade ()2U|) and (|19j) . We stress that in Theorem El we use a different assumption from 
USC which is typically weaker than the latter. Moreover, we use Lyapunov SP-UAS which 
is typically stronger than SP-UAS. In particular, in Theorem El we may have that the right 
hand side of (|2T)|) is discontinuous, which is in general excluded from Theorem El because of the 
USC condition. Note also that since appropriate converse Lyapunov theorem does not exist for 
parameterized systems (see Remark |2J), we only state sufficiency results in Theorem El Finally, 
we note that the proof of Theorem El is inspired by the proof of [2H1 Lemma 2] (see Section EJ). 

Theorem 2 Suppose that the solutions of the system Alfy) with input z are USC (UGC). Then, 
the system < fT%j) . < TH^) is SP-UAS (UGAS) if and only if the following conditions hold: 
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1. The system IJ® is SP-UAS (UGAS); 

2. The system is SP-UAS (UGAS); 

3. The system hlty . M$) satisfies the property USB (property UGB). 

□ 

In order to state our second main result we need the following assumption. 

Assumption 1 There exist 72 G N ' , 71,73 G /Coo o,nd T* > such that for all £ G M. n , k > 
and T G (0, T*) we nave: 

\f T (k,x,z)\ < 71(1^1) 
|/t(A;,z,*)-/t(M,0)| < T l2 (\x\) l3 (\z\) . 

□ 

Theorem 3 Suppose that fx of the system \lty satisfies Assumption^ Then, the system \lty) . 
hlfy) is SP-UAS (UGAS) if the following conditions hold: 

1. The system Wty is Lyapunov SP-UAS (Lyapunov UGAS); 

2. The system \J$) is SP-UAS (UGAS); 

3. The system J7%J) . \~T9\) satisfies the property USB (UGB). 

□ 

For the sake of clarity we present the proofs of the main results in Section We provide 
only the proofs for the more general case of SP-UAS since the global versions follow by removing 
the restriction on the size of the domain of attraction and restricting the neighborhood of the 
origin to the origin itself (i.e., considering v = 0). 

Remark 3 While the proof of Theorem^ is more direct since it does not require the existence of 
a Lyapunov function, Theorem^ is very important since the existence of a Lyapunov function 
is often helpful for controller design. Indeed, a Lyapunov function may allow us to improve 
the transients of the sampled- data system by redesigning a continuous-time controller in an 
appropriate way and using our results. 

Two interesting corollaries for non-parameterized systems follow directly from our results. 
These results are interesting in cases when the exact discrete-time model of the plant can 
be computed and we do not have to appeal to Theorem ^ The appropriate definitions and 
assumptions for non-parameterized systems, 

x(k + l) = f(k,x{k),y(k)) (28) 
y(k + l) = g(k,y(k)), (29) 

are obtained easily from global definitions for parameterized systems by setting T = 1 and they 
are not repeated for space reasons. 
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Corollary 1 Suppose that f is continuous in x,y, uniformly in k > 0. Then, the system 
is UGAS if and only if: 

1. The following system is UGAS 

x (k + 1) = f(k,x(k),0) ; (30) 

2. The system is UGAS; and 

3. The system is UGB. 

□ 

Corollary 2 Suppose that Assumption^ holds for f in fg£j) . Then, the system $EB), is 
UGAS if: 

1. The system A3U\) is Lyapunov UGAS; 

2. The system is UGAS; and 

3. The system » UGB. 

□ 



We stress that UGB is in general difficult to check. In jzo] we present several sufficient 
conditions for this property to hold and which are inspired from EE]- For the sake of 
completeness, we close this section with a result which is representative of the type of conditions 
given in those references and which we will use later on in the case study. To that end we 
introduce first the following technical hypothesis. 

Assumption 2 Suppose that there exist ai,a 2 ,(p € /Coo, 7i? 72 £ A/", c, T* > and for each 
T e (0,T*) there exists V T : M> x W 1 * -> M> such that for all x eR"v G IR™ 2 , k > and 
T E (0,T*) we have that : 

&i(\x\) < V T {k,x) < a 2 {\x\) + c (31) 
V T (k + l,f T (k,x,z))-V T (k+l,f T (k,x,0)) < T^(\z\) V (V T (k,x)) + TU\z\) (32) 
V(k + lJ T {k,x,0))-V T {k,x) < (33) 

ds 



i <p(s) 



oo . (34) 



Proposition 2 Consider the system hlty) with input z and under Assumption® If furthermore 
the solutions of the system hlty satisfy the summability condition 

Tj2^(\Mk,k ,z )\)<p(\x \), (35) 
with some p G /Coo and p(s) := 7i(s) + then, the system is UGB. □ 
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The proof follows along the lines of |2~3*1 Lemmas 1 and 2] and therefore is omitted here 1 . 

The proposition above establishes some interesting links with conditions used in the litera- 
ture in the context of continuous-time systems to prove UGB. The condition (|34|) restricts the 
growth of (p. In particular, it holds when ip(s) = s; this situation was considered for instance in 
[3*H] with 7i(s) = 0. Earlier results using similar conditions are found in [2B] and in the context 
of forward completeness already in (HH] • 

In the particular case when 7 1; 72 G /Coo an d c = 0, one can prove that there exist /C^ 
functions a\, «2 such that the systems trajectories satisfy 

fc-i 

ai (\(f) x T (k,k ,x ,u;( ko , k) )\) < a 2 (\x \) + Tj2KH ,k)\), Vk>k o >0. (36) 

The property of the system by which it satisfies (|3l)|) is called integral input-to-state neutral 
stability (ilSNS) and has been introduced for non-parameterized discrete-time systems in 
See also [32]; for the original definitions of integral input-to-state stability in the context of 
continuous-time systems. 

The summability condition (J33j) which imposes a minimal convergence rate for the system 
(fTTJj) is tight. However, we are not aware of any proof of necessity. This condition is in general 
hard to check since it is trajectory-dependent. As far as we know such condition was firstly 
used in the context of continuous-time cascades, in [SZj. Sufficient conditions for this property 
to hold are also given in [2H] along the lines of the results presented in [2] . 

Finally, we shall mention that the condition (|3*2*|) holds for instance (under suitable growth 
conditions on Vr{k, •)) if GT{k,x, z) has linear growth in x for each fixed k and z. This type 
of conditions have been extensively used in the continuous-time context to avoid the so called 
peaking phenomenon. See [1U] and references therein. 

4 A case study: tracking control of the unicycle 

In this section we revisit the problem of tracking control of a mobile robot of the unicycle type. 
This problem has been thoroughly studied in the continuous-time context via many different 
approaches (see ^Hj for a survey; for a more recent text with an updated list of references see 
(2D])- To illustrate the utility of our results we will revisit the cascades approach used in [HE] 
for a 3 degrees of freedom cart. The results may be extended to higher dimension systems, 
following for instance |2*U] . 

While the problem setting is the same as considered in the continuous-time context, we will 
see that the proof techniques employed in the discrete-time case are quite different. For instance, 
since we deal with approximate discrete-time models, some important structural characteristics 
are lost. Hence, we believe that the proofs of this section are interesting in their own right. 

According to ^3] the context of the problem can be set as follows. We have a mobile robot 
with two directional wheels and two "fixed" wheels and whose motion is described by 

x — v cos 9; y = v sin 9\ 9 = uj , (37) 

where x, y are the Cartesian coordinates of the center of the axis joining the directional wheels 
and 9 is the orientation angle of the directional wheels. The robot is required to follow a 
trajectory generated by an exosystem, i.e., a fictitious "reference robot" with kinematics 

x r = v r (t) cos# r ; y r = v r (t) sin^; 9 r = u r {t) (38) 
1 For review purposes we provide it in the Appendix. 
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where v r (t) and u r (t) are given reference velocities. Then, the tracking errors satisfy the set of 
equations (see O Lemma 1]) 

x e = u>y e — v + v r (t) cos9 e ; y e = —ux e + v r (t) sin# e ; 9 e = u> r (t) — ujv , (39) 

where (-) e := (-) r — (•). The system is velocity-controlled, i.e., the control problem reduces to 
finding control inputs uj and v (which also correspond to the actual angular and linear velocities 
of the cart) such that the origin of is UGAS. 

There are numerous solutions to this problem in the context of continuous-time (e.g. [TT|l2l1] 
and j20j for a recent literature review). Here, we will revisit the cascaded approach proposed 
in whose main feature is that the control laws are linear. To illustrate and motivate our 
results we will proceed to solve the same problem with a linear time-varying discrete-time 
controller which we will redesign based on the Euler-discretization of the error dynamics, 

x e {k + 1) = x e (k) + T[ouy e (k) — v + v r {k) cos9 e (k) } 

y e (k + l) = y e (k) + T[-ux e {k)+v r {k)sm9 e {k)} (40) 

9 e (k+l) = 9 e (k) + T[u r {k) -co}. 

Thus, our control problem consists of designing v and uj such that (}40|) is UGAS. 

We will solve this control problem following a similar approach to that of where it was 
shown using results for continuous-time cascaded systems, that the system (JBTJj) in closed loop 
with v = v r (t) + a2X e and uj = uj r (t) + a\Q e is UGAS for appropriately chosen a± and a^- 

The controller structure that we use 



uj := u) r {k) + ai9 e (k); v := v r {k) + a2X e {k) + T-& (41) 

where a±, 02, u r (k) and v r {k) come from the continuous-time control law proposed in and 
d is an extra control input which depends on fc, i e and y e and that we will design with the aim 
of improving the system's performance. More specifically, the motivation for the control laws 
above is that as in the continuous-time context, the closed loop system 

x e (k+l) = (1 - Ta 2 )x e {k) + Tu r {k)y e {k) - T 2 ti 
y e (k + l) = y e (k) -Tuj r (k)x e (k) 

V v ' 

F 1T (k,x(k)) 

V V (42) 
T[ a\9 e {k)y e (k) — v r (k) + v r (k) cos9 e (k) ] 

T[ —ai9 e {k)x e {k) + v r {k) sin 9 e (k) ] 



G T (k,x{k),z{k)) 

9 e (k + l) = (1-T ai )9 e (k) =: F 2T (k J x(k)), (43) 

where z := 9 and x := col[x e , y e ], has a cascaded structure. 

Hence, the control laws (}4T|) are designed with two main ideas in mind: 1) to have as simple 
as possible controllers; 2) that the closed loop system verifies the conditions of our main results 
for cascades. More specifically, notice that the bottom subsystem (j4^|l is independent of x e and 
y e and is UGES for values of a\ sufficiently small (T*a x < 1). Also, the interconnection term 
Gr(k,x,z) is linear in x e and y e . Hence, our results suggest that we only need to design $ as 
a function of x e and y e only so that the zero-input (i.e., with Gt = 0) subsystem in ()42j) be 
UGAS (or possibly, UGES). 

Notice that in the particular case that $ = we obtain the emulated (discretized) continuous- 
time control law. However, as we will illustrate below, when carefully defined this extra degree 



13 



of freedom in the control design allows to improve performance and, on occasions, to enlarge 
the domain of attraction with respect to that of the emulated continuous-time control law. 
Thus, our control scheme can be regarded as a redesign of the cascaded-based continuous-time 
controller of [SB] • Simulation results at the end of this section will illustrate the advantages of 
this approach. 

Here, we will establish UGES of the zero-input system based on a property of persistency 
of excitation. However, we will need a specific reformulation of this property 2 within the 
framework of discrete-time parameterized systems. This is introduced next. To compact the 
notation, in the sequel we will use u Tk := ui r {k). 

Definition 5 (PE) Let u r : Z>o — > Z be a junction produced by sampling a function ip : 
M>o — > M at rate T. The function uo r is said to be persistently exciting (PE) if there exist 
positive numbers fi, L and T* such that for all T £ (0, T*) and all j > 0, 

r£<>A*. (44) 

k=j 

We are now ready to present our main result of this section. 

Proposition 3 Consider the system fi^ffl in closed loop with \4l\j . Assume that 

1. There exists T,wm > such that for all k > and T £ (0, T) 

\Vr h \ , \uJr k \ , J Tf, 1 } < W M ■ (45) 

2. The signal u rk is PE. 

Then, there exists a 2 > such that for all K,a% > and i9(fc,x) with \$(k,x)\ < K\x\, the 
system is UGAS. 

Proof. It follows by invoking Theorem El Firstly, we see that Assumption Q holds trivially 
in view of item 1 of the proposition. To see more clearly, notice that (|45J) implies that fr{-, ■, •) 
is continuous and uniformly bounded in the first argument. Item 1 of the proposition also 
implies that there exists c > independent of T, such that for all T £ (0,T*), we have that 
\G T (k,x,z)\ < Tc\z\ (\x\ + 1). 

Secondly, it is evident that the origin of (f4l?J) is uniformly globally exponentially stable 
for any a\ and any T £ (0,T*) where T* > is such that 1 > a{T* > and therefore the 
trajectories </>y(-, •, •) are uniformly summable. 

It is left to prove that the unperturbed dynamics x(k + 1) = Fi T (k,x(k)) in (J4"2)) is UGAS 
and that (J4~2|) is UGB. As a matter of fact we will show that the zero-input system in (}4~2|) is 
Lyapunov UGES. 

Proof of UGES of x{k + 1) = Frr(k, x(k)): Consider the function Vr(k, x) := \x\ 2 — eu rk l x e y e 
with e := a y + T and a y > 0. Observe that this function is positive definite and radially 
unbounded for sufficiently small a y , T* and wm] indeed, we have that 

ci |x| 2 < Vr(k, x) < C2 \x\ 2 (46) 
2 With respect to the original one from 0. 
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with c\ := (1 — 0.5(a y + T*)wm) and c 2 := (1 + 0.5(0^ + T*)w M ) which are clearly independent 
of T. We assume that a and T* are sufficiently small so that ci > 0. Next, we compute: 



AV T := V T {k,F 1T (k,x))-V T (k,x) 

-T{2a 2 -eu 2 rk )x 2 e -eTul k y 2 e 



+T 2 x e ([a 2 , + u 2 k - ea 2 ]x e - [2a 2 u rk - eu* h ]y e - [2(1 - a 2 T) + eu 2 r T}§) 
+T 2 y e (uj 2 k y e + eu) rk d) + T 4 $ 2 + ex e y e (u rk - LO rh _ x + uj r . k a 2 T) . (47) 

Under the assumptions of the proposition, ex e y e (u rk — uj rk _ 1 + u rk a 2 T) < (1/2) \T 2 y 2 + 
e 2 u;| / (l + a 2 ) 2 x 2 ]. Define a x := a 2 - (a y + T*)w 2 M - (l/2)e 2 w 2 M (l + a 2 ) 2 which is positive for 
sufficiently small values of e and sufficiently large values of a 2 . Also, since \'d(k,x)\ < K \x\ 
there exists K± > such that 

^<-Kx 2 + W 2 fe y 2 )+TAMx| 2 , (48) 

for all T G (0, T*), x G 1R 2 and fc > 0. 

Let us introduce the following auxiliary function: 

oo 

W T (k,x) := -Tj2e (k ~ l)T ^y 2 e 

i=k 

for which we claim the following (for the proof, see the Appendix). 

Claim 1 Suppose that the signal ou rk is PE. Suppose also that there exists wm > such that 
for all i > we have \u n \ < Wm- Then, there exist strictly positive numbers T*,C3,C4,K 2 ,a y 
such that for all T G (0, T*), k > and z e G M. 2 we have 

-c 3 y 2 e < W T (k,x) < -c iV 2 e (49) 

AW T 9 9 9 9 , x 

-j^ < u 2 k y 2 e - a y y 2 e + A 2 a; 2 . (50) 

□ 

Let now T* be generated by the claim above and such that (|4T)|) and ()48j) hold. Then, we can 
complete the proof by showing that there exists e > and T > such that Uxik, x) := Vr{k, x) + 

eWr(k,x) is a Lyapunov function that proves UGES. Indeed, let e := min a^j and 

T := min |t*, min ed y } |, where T* > comes from the proof of the proposition and 
ct y comes from (|5T)|) . Then it is easy to show that Ut satisfies 

L 'x| 2 < Uxik.x) < c 2 \x\ 2 (51) 



2 



AU T 

< — It 

T 



< -cz\x\ , (52) 



where C3 = | min{^-, ea y }. This completes the proof invoking standard Lyapunov arguments. 

Proof of UGB : We invoke Proposition It is worth recalling to avoid confusion in the nota- 
tion, that the state x in Proposition El corresponds here to x = col[x e , y e ] and the input z in 
Proposition |21 corresponds here to 9 e . Hence, we proceed to verify the conditions of the proposi- 
tion with f T (k,x,z) := Fi T (k,x) + GV(/c, x, z) as defined in (|4^|) and Vx(k,x) = Uxik.x). The 
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bounds flHU) and ® hold from (fET jl and The conditions ((33) and flHl hold with p(s) = s, 
72 = c? 2 s, and 7i(s) := d 3 s, d 2 , d 3 > 0. This is because Uxik, ■) is quadratic and Gx(k,x,z) 
contains terms of linear growth in x for each fixed k and z and trigonometric functions, which 
can be over-bounded by a linear function of x. Also, GT{k,x,-) can be over-bounded by a 
linear function for each fixed k and x. Finally, holds with a linear function p(s) := p s 
since p(s) is in this case a non decreasing function of linear growth and (j)^(k) decays uniformly 
exponentially to zero. ■ 

Remark 4 It is worth remarking that in the proof above we could also have invoked Theorem^ 
Notice that to show that the system Alty) is UGC with input 9 e we may appeal to Lemma[I\ob- 
serving that the system x(k + l) = fT(k,x(k),9 e (k)) with f T (k,x,9 e ) := Fi T (k, x) +Gix(k, x, 6 e ) 
is linear in the state x = col[x e , y e ] and for each A > 0, we have that fT(k,x,8 e ) is globally 
Lipschitz in 9 e uniformly for all k > and x such that \x\ < A, with a Lipschitz constant of 
the form L = L'T with V depending only on A. 




Figure 1: Tracking errors for x and y and the control input v given by (j41)l . (|5Hj) . 

Now we illustrate how we can use Proposition El to improve the performance of the system 
with the redesigned controller. Since the correction d can be chosen arbitrarily, we can choose 
it so that negativity of AV in ()47|) is enhanced. Indeed, a closer inspection of the difference 
equation (JTTj) shows that one such choice is 

vik.x) :— — ; — — 53) 

V 7 2(1 -a 2 T) +eu* k T V ; 

with e = a y + T. 

The simulations that we present next show that the performance is considerably improved. 
We have simulated the system above in Simulink™ of Matlab™ with a 2 = 70, a± = 10, 
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w r (k) = 20 sin(kT), T = 0.01 and a y = 2 — T. We show the results in Figure ^ We show only 
the responses for the states x e and y e as well as v since these are the only variables affected 
by the additional input T§. We show simulations for the system's response with d = 0, with 
$ = 0.5 [ (al + — ea<i)x e — (2a2^ rfc — £UJ f h )y e ] and for $ as defined in (IH3j) . The best apparent 
performance is for the latter. 

It is also clear from the plots, that even though the correction d is linear in the state and 
actually (d ~ 0(1) \x\), this correction is not comparable to "adding gain" to the control input. 
Notice that in this case the resulting control effort is actually smaller than in the case of the 
continuous-time based controller (i.e., when # = 0). 

5 Proofs of the main theorems 

All results are only proved for semiglobal-practical stability properties since global results follow 
the same steps with minor changes. 

5.1 Proof of Theorem [2] 

Lemma 1 Suppose that the solutions of the system hltyl with input z are USC (UGC) and the 
system is SP-UAS (UGAS), with the function (3 X G ICC Then, for any pair of strictly 
positive real numbers (A, v) there exists T* > (there exists T*) such that for any rj G (0, A] 
(any rj > 0) and any strictly positive real numbers (e, L) there exists fi > such that for all 
inputs z(-) with \u z \ < \x, k Q > 0, x(k ) = x Q with \x \ < f] and T G (0, T*) the following holds: 

\(j>T(k,ko,x ,uf koik) )\ <max{/3 x (|x |,(fc-fc )T),i/} + e, (54) 
for all k G [k a , k Q + £l,t\- D 

Proof. For any (A, v) there exists T* > such that for any 77 G (0, A] and any strictly 
positive real numbers (e, L) there exists fi > such that for all inputs z(-) with \u z \ < fi, all 
k a > 0, x(k ) = x with |x D | < rj and T G (0,T*) the following holds: 

\4>r{k, k ,x ,uf kojk) )\ < \^(k,k o ,x o ,0)\ + \^(k,k ,x ,ujf ko>k) ) - (j)^(k,k o ,x o ,0)\ 

Hence, (|B^) follows from (J23j) and the KL bound on \4>^(k, k a , x a , 0)|. 

Lemma 2 Suppose that all the conditions of Lemma Q hold. Then, for any strictly positive 
(A, v) there exists T* > such that for any e > and r\ G (0, A] there exists fi > such that 
for any input z(-) with \u z \ < fi, k > 0, x(k Q ) = x with \x Q \ < A and T G (0,T*) we have 
that the solutions of the system with input z satisfy the inequality \54\j for all k > k Q . □ 

Proof. Let (3 X G ICC come from SP-UAS stability of l(2njl . Let (A, v) be given. Without loss 
of generality assume that A > v and 

AzO,0)>s, Vs>0. (55) 

Let v G (0, 1) be such that 

P x (v,0) + v/2< v . (56) 
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Note that (|55|) implies that v < v. Let (A, v/2) generate T* > via Lemma E] and assume 
without loss of generality that T* < 1. Let e > and r\ G (0, A] be given. Let b\ > be such 
that 

P x (6 u 0) + 6i/2<e. (57) 

Let L > 1 be such that 

/^(max^A, < y ■ (58) 

Let ei := min{i>/2, Si/2} and let e± and L generate /x > via Lemma^ In the rest of the proof 
we consider arbitrary fixed k a > 0, x(k a ) = x a with \x Q \ < 77, \u z \ < fi and T G (0, T*). 
In order to simplify the notation, for i G Z>i we introduce: 

fcj := k + i£ LtT 

Xi := (j) T (ki,ko,x ,U{ ko ^) . 

Time interval [A; ,fci]: From Lemma Q (JEEJ), (jHHj) we can write 

|0y(/c, fc ,Zo,^f fco)fe ))| < max{/3 :c (|x | , (k - k )T),u/2} + e x 

< max{f3 x (\x \,(k-k )T),is} + e (59) 

for all k G [fc , fci]. Moreover, note that 

T < 1 L — 1 < T£ LtT . (60) 

Using, (JHHJ), ij^Hjl and our choice of e\ we can write 

\ Xl \ = \cj)*(ki,k ,x ,uf koM) )\ < max{/3 x .(|x | , T£ L>T ), v/2} + ei 

< max{/? x (?7, L - l),i>/2} + ei 

< max{5!/2,i>/2} + min{(5i/2,i>/2} 

< max{5i, u} . (61) 

Time interval [fci,/^]: Since |a?i| < max{5i,£}, we consider two cases. 
Case 1: \x\\ < v. In this case, using Lemma [H (J5T)j) and (|55j) we can write 

\^(k,k ,x ,u( ko)k) )\ = \^(k,ki,xi,u( kuk) )\ 

< maxj^dxil ,0), z>/2} + ei 

< max{/3 x (i>,0),£/2} + min{i>/2,5 1 /2} (62) 

< p x (y,0) + i>/2 < v , 

for all k G [fci, fc 2 ]. Moreover, using (|6T|). the fact that u < 1, (j3HJ), dS0|) an d our choice of e x we 
can write 

|0r(^2,&o,a; o ,c^ o;fc2) )| < max{/3 x (|xi| ,T£ i>T ),i>/2} + e x 

< max{/5 x (P, L — 1), z>/2} + ei 

< max{/3 a .(l,L - 1), z>/2} + e 1 (63) 

< max{<5i/2, i>/2} + min{5i/2, i>/2} < max{5i, P} . 

Case 2: \x±\ < Si. Using Lemma Hand (jHTjl we can write 

\(f)^(k,k ,x ,ujf kotk) )\ = \<f)^{k,ki,xi,u^ k) )\ 

< max{/3 x (5 1 ,0),z>/2} + min{z>/2,(5 1 /2} (64) 

< max{/3 x (<Ji, 0) + Si/2, u} < max{e, u} , (65) 
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for all k G [ki, k 2 }. Also, in a similar manner as before it follows that 

\(j)%(k 2 ,k ,x ,L}f koM) )\ = \</%(k 2 , fci,xi,wf fclifc2) )| 

< max{p x (5 1 ,T£ LtT ),i>/2} + e 1 (66) 

< max{&(<$i,L-l),i//2} + ei (67) 

< max{<5i/2,i>} + min{i>/2,5i/2} 

< max{5i, u} . (68) 

Time intervals [ki, fcj+i], z > 1: Using similar calculations it can be shown by induction that 
for all integers % > 1, 

\(j)^{k,k Q ,x ,uf koik) )\ < max{e, u} , Mk G [ki, k i+1 ] (69) 
|0r(/ci,/co,^o,wffc o , fcl ))| ^ max{5i,z>}. (70) 
The proof follows from ()59|) and ()69)) by noting that 

\4>^(k,k ,x ,u^ ko ^)\ < max{e, u} < e + u 

< m&x{j3 x (\x \ , {k — k Q )T), u} + e , VA; > ki . (71) 

Lemma 3 Suppose that all conditions of Lemma[I\ hold and (3 X G ICC comes from SP-UAS 
(UGAS) of I120\) . Then, there exist c x , c u > 0, 7 G /C and for any v > there exists T* > swc/i 
taai for all k Q > 0, x(k Q ) = x with \x \ < c x , \u z \ < c u and T G (0, T*) the following holds: 

\(/)^(k,k ,x ,ujf ka>k) )\ < max{/5 x (|x | , (k - k )T), u} + 7 (|^f feoiA:) |) (72) 

for allk>k o >0. □ 



Proof. The proof of this fact follows closely the proof of fSl Lemma 3.3]. Since all conditions 
of Lemma d hold, the conclusion of Lemma El holds. Let [c x ,v) generate T* > via Lemma 
[0 Let 7] = c x and for a fixed e > 0, let /2(e) be the supremum over all applicable /i's. Then, 

< /i(e) implies that (f^H) holds for all k > k Q > and if fii > Ji(e) then there exists 



to 



[feo.fe) 



^ffco.fc) 



< /ii and 



0f,(/c, fc , i , oj? ka k ^) > max{f3 x (\x \ , (k — k )T), u} + e. /i(e) is positive and 
nondecreasing but it is not necessarily continuous. Chose ( G /C such that C(r) < kfi(r), with 
/c G (0, 1). Let 7(7") = C~ 1 ( r ) an d note that 7 G /C. Let c u := lim^oo C( r )- 

Let |x | <c x ,Te {0,T* 
have that (72) holds for all k > k >0. 



Given uf ka k ^ with 



a; 



[fco.fc) 



< c u let e = 7 



Then, we 



Remark 5 Tne following fact was proved in \3$ : for any j3 G ICC and any c > there exists 
$ G /CX such that 

P(s,t) <(3(s,t + c) , Vs,te [0,oo). 

27m can fre further strengthened in the following manner. It was shown in \4'Jj that given any 
(3 G /CX, there exist a G /Coo, K G K, such that /3(r,t) < cr(fi:(r)e _ *) Vr, t > 0. Consequently, 
given any (3 G /CX and a nondecreasing function L(-), there exists (3 G /CX swca iaat 

P(s,t)<p{s,t), Vie [0,L(s)] . 

To show this, let (3 generate a, k G /C^ as above. Then (3(s,t) := a(K(s)e L ^e~ t ) proves the 
claim. □ 
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Lemma 4 Let £ = {x T z T ) T . Suppose that there exist /3i,/3 2 ,/3 3 G ICC, 7 G /C and c > snc/i 
that for all u x > 0, tnere exists T* > smc/i that for all max{|n(/c )| , |f(fc )|} < c , T G (0,T*) 
and k > k Q > we have 



\x(k)\ < max{p 1 (\x(k )\ , (k - k )T), v{\ + 7 ( sup |n(t)| 
|u(fc)| < max{/3 2 (|e(A; ) 1,(^-^)7),^} , 

\z{k)\ < m^{^{\z{k Q )\,{k-k )T),u y } , (73) 

where £ := (x T z T ) T . Then, there exist (3 G /(X and ci > suc/i that for any v > tnere exists 
T* > snca taat /or a// f (fc ) = £0 wita |£(fc )| < ci, T 6 (0, T*) and k > k Q > we have 

\C(k)\<msix{P(\Co\,(k-k )T),u} . 

In particular, we can take: 

/3(s,t) = 4Af2/3 1 ( S ,t/2) + 2 7 (/3 2 ( S ,0)),t/2 N )+47(/3 2 (s,t/2)) + 2/33(s,t) , (74) 



where all are generated via and an arbitrary fixed c > using Remark^ Moreover, if all 
the inequalities hold globally, then the conclusion holds globally and we can take: 

(3(s,t) = (&(s,i/2)+ 7 (&(s, 0)), t/2) + 7 (&(s, */2)) + fo(s, t) . (75) 

□ 

Proof. Let c > be fixed and let $ and c generate via Remark Let c come from the 
Lemma and let v be given. Let ci, z/i be strictly positive numbers such that 

c > max{2/3 1 (c 1 ,0) + 2 7 (/3 1 (c 1 ,0)),/3 2 ( Cl ,0),c 1 } 

f > 7o(^i) , (76) 

where 70 (s) := 4/3i(2s + 27(5), 0) + 47(5) + 2s. Let i^i generate T* > via the Lemma. Let 
T* = min{T*,c}. Define a piecewise constant function x(-) as x(t) = x(k),t G [kT, (k + 1)T); 
z(-), u(-), £(•) are defined in the same manner and can be regarded as functions that depend on 
continuous time t e 1. From the definition of $ and our choice of T* we have @i(s, t + T*) > 
Pi(s,t + c) > Pi(s,t), Vs,t- 

Let |f (t ) I < ci and T G (0,T*). Then the following holds: 

|x(t)| < max{^i(|x(t )| , t - t ), v{\ + 7 sup \u(r) 

\re[t ,t] 

\u(t)\ < max{&(|£(t )|,t-toW , 

\z(t)\ < max{^ 3 (k(^)|,t-to),^i} , (77) 

for all t > t > 0. Using the same technique as in fSJ pg. 221-222], we first find a bound for 
x(t) starting at ^y 2 - and then a bound for x (^y^) starting at t . Using the same procedure as 
in ^Hl pg- 221-222], these inequalities together with the following facts 

max{a, b} + max{c, d} < max{2a + 2c, 26 + 2d} , Va, 6, c, d > 
sup |u(r)| < /3 2 (|f(t o )|,0) 

re[to,(t+to)/2] 

sup |n(r)| < /3 2 (|e(t )|,(t-t )/2) 

re[(*+*o)/2,t] 
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generate the inequality 

\C{t)\ <max{/3(|e(to)|,t-to),7o(^i)}- 
This inequality, together with our choice of v\ and the fact that £(kT) = implies that 

\S{h)\<max{/3{\S{k )\,{k-k )T),v} 

which completes the proof. 

Proof of Theorem [2J Let j3 x G ICC come from item 1 of the Theorem. Let (3 Z G K.C come from 
item 2 of the Theorem. Let k G /Coo and c > come from item 3 of the Theorem. Let 7 G /C 
and c x ,c u come from Lemma El Let c z > be such that (3 z {c z , 0) < c„ and := mm{c z ,c x }. 
Let /9 G /C£ and c\ > be generated using /?i = /?2 = /?3 = /3 2 , 7 and Co = via Lemma 

HI Let «i(s) := + lj k(s). Let Lj : 1R> — » IR> , z = 1,2 be continuous nondecreasing 

functions with Lj(0) = 1 such that for all A > we have 

P Z {^L X {&)-1) < Cl (78) 
/^(A^oz^A) -1) < ci/2. (79) 

Let L(s) := Li(s) + L 2 o /«i(s). Let /3 and L(s) generate j3 via RemarkEl Finally, we define 

P(s,t) :=max{^(K 1 (s),t),^(s,t),fi:i(s)e Ll(s)+i20Kl(s) e^} . 

Let (A, v) be given. Let v\ := min{ci/2, z/}. Let (k\(A), V\) generate T* > via Lemma ^ let 
(A, ui) generate T 2 * > via the item 2 of the Theorem; let v\ generate T 3 * > via Lemma El let 
v x > generate T 4 * > via Lemma 13 Let T* : = mm{T* , T 2 * , T 3 * , T 4 * , 1 } and T G (0,T*). Let 
> 0, £(k ) = £ Q with |£ Q | < A. We consider two cases. 

Case 1: If |£ Q | < c±, then by a direct application of Lemma El and |U we have that 

4(Mo,£o)| < max{P{\£ \,{h- k^T),^} < max{/3(|£ | , (k - k Q )T), u}, (80) 
for all k > k Q > 0. 

Case 2: If |£ Q | G [ci, A], then since |£ Q | > c 1; we can write using item 3 of the Theorem: 

4(Mo,£o) < «(ie |) + c < «(|6I) + -7^ T «(|6I) = « 1 (|6I) , (81) 

for all fc > k D > 0. Hence, we can write that for all k G [k D , k Q + £l-l,t + £l 2 ,t]' 

4(Mo,£o)| < ^(led) < ^(leoDe^d^lH^^idfoDe-ffc-feJT . (82) 

Using definitions of £l u t an d the fact that T < T* < 1, we obtain that 

- 1 < T^, T < 2 = 1,2. (83) 

Moreover, the definitions of Li, v we get that for all k > k + ^Li,T) 

|<^(A;,A; ,z )| < max{&(A, TY LljT ), u} < max{/5,(A, L : (A) - 1), v} < c x . (84) 

Furthermore, by letting k\ := k Q + £ LlsT , k 2 := k Q + 4 1)T + £ L2iT , x\ := ^(h, ko,x ,u^ koM) ) 
and x 2 = 4>x(k2i k , x , ujf ka k ^) we obtain that 

\<p^{k 2 ,k ,x ,u;f koM) )\ = \<^(k 2 ,k 1 ,x u uf klM) )\ 

< vaax.{/3 x (\xi\ i (k 2 -k 1 )T),Vi} + e 

< max^/^A), T^ 2)T ), i/J + Cl /2 

< max{/3 :c (re 1 (A), L 2 o « X (A) - 1), 1^} + c x /2 

< max{ci/2,ci/2} + ci/2 = c x . (85) 
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Hence, using f)84|) . ((53)) , Lemmas H3 and H] and the construction of /3 we have that for all k > k 2 



(j)j^^k, k , £ 



4> T (k, k 2 ,& 



< max{0(|£ 2 | , (A; - fc 2 )7Vi} 

< max{^(«i(|^ |), (k-k - £ Ll)T - £l 2 ,t)T), vi} 

< max{^(Ki(|^ |), (k - k Q )T -L x - L 2 ), i>i} 

< max{^(Ki(|eo|),(A;-A;o)r),i/} , 

which proves the result with the defined (5. 



5.2 Proof of Theorem [3] 

We again prove only semiglobal results since global results follow the same steps with minor 
changes. 

Lemma 5 Suppose that all conditions of Theorem^ are satisfied. Let Vt come from item 1 of 
the Theorem and define 



y(k, := V T ( k + 1, f T (k, x,z))- V T { k + 1, f T (k, x, 0) ) 



(86) 



Then, there exists ft y G K.C such that for all (A, v) there exists T* > such that for all k a > 0, 
£(k ) = £o with |£ Q | < A and T G (0,T*) we /iave t/jai 



(f) T (k, fc , £,)) < T ■ max{^(|^ | , (fc - fc )T), ^} . 



/or a// k > k 



(87) 
□ 



Proof. Let (A, v) be given. Let 71,73 G /Coo, 72 £ A/" and T* > come from Assumption [TJ 
Let k and c come from item 3 of the Theorem and let A generate T 2 * > via the same item. 
Let L G M and A x := 7 1 (k(A) + c) generate T 3 * > via item 1 of the Theorem. Let (3 Z G ICC 
come from item 2 of the Theorem. Define 



/3 y {s,t) 



L o 7i(ac(s) + c) • 72(/c(s) + c) • 73 o ^(s, t) 



7 3 



Io 7l ( K (A)+c)- 72 WA) + c) 



Let (A, fi) generate T 4 * > via the item 2 of Theorem. Let T* := minjT^, T 2 *, T 3 *, T4 } and 
consider arbitrary fixed k Q > 0, £(fc ) = £0 with |£ Q | < A and T G (0, T*). Then, using item the 
definition of y, Lipschitz condition on Vt and Assumption ^ we can write 



\y(k,£)\ = \V T {k+l,f T {k,x,z))-V T {k+l,f T {k,x,0))\ 

< T ■ L (max{/ r (fc, x, z), f T (k, x, 0)}) \ f T (k, x, z) - f T (k, x, 

< T.Lo 7l (|e|). 7 S > (W)-7sW , 



(88) 



for all |£| < Ai, T G (0, T*). Finally, using items 2 and 3 of the Theorem and the definitions of 
/3 y and z/i we can write 

y (k,(f) T (k,k ,£ )\ < T-Lo 7l ( <f) T (k,k ,£ ) ) ■ 72 ( fc , £0) J ■ j 3 (\<f) T (k,k ,z )\) 

< T-Lo 7i (k(|£ |) + c) • 72(«(|e o |) + c) • 7 3(max{/5 z (|^| , (A; - k )T), v,}) 

< T-m&x{(3 y (m,(k-k )T),v} , (89) 
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for all k > k . 

Proof of Theorem We only show that there exists (3 X G JC£ and 7 G K, such that for any 
(A, v) there exists T* > such that for all k Q > 0, f (fc Q ) = £ Q with |f | < A and T G (0, T*) we 
have that: 

|0^(/c,/c o ,x o ,u; 2 )| <max{^ x (|xo|,(fc-fc )T),i/} + 7(|c;f fco)fc) |) , (90) 

for all k > fc , where m(/c) := y ^ k ^T(^ k °^°)) an( j w z j g a se q Uence Q f u(k)'s on the interval [fc , fc). 
Then, the conclusion of the Theorem follows directly from Lemma item 2 of the Theorem 
and Lemma El 

Let ai, 0(2,0:3 come from item 1 of the Theorem. Define «(•) := -03 o o 2 ^ 1 (-) and consider 
the differential equation: 

C = -a(v) ,C(0) = Co ■ 

Without loss of generality we can assume that a is locally a Lipschitz function (see fSl footnote 
on pg. 139]) and from j3H Lemma 2.5] we have that there exists (3 G /C£ such that for all < ( 
the solution of the differential equation exists, it is unique and satisfies: 

C(t)=p(Co,t-t ) , vt>t o >0. 

Define a(s) := a _1 (s) + s and: 

P v ( 8 ,t) := a^{(](a 2 {s),t)) 

7 (s) : = o" 1 oa(2s) . (91) 

Let (A, z/) be given. Let Ax := k(A) + c, where k, c come from the item 3 of Theorem. Let 
v\ := o Oi(v). Let (Ai, ^1) generate T* G (0, 1) be so that the items 1 and 3 of Theorem 
hold. Consider arbitrary fixed k Q > 0, £(k ) = £0 with |£ | < A and T G (0, T*). Note that 
the item 3 of Theorem guarantees that the solutions of ()18j) . (j!9j) exist for all fc > A; > and 



satisfy 



<Ar(/c, fco> £0) 



< Ai. We can write 



V T {k + l,f T {k,x,z))-V T {k,x) = V T {k + l,f T {k,x,0))-V T {k,x) + 

V T ( k + 1, f T (K x,z))- V T (k + 1, f T (k, x, 0) ) 
= V T (k + l,fT(k,x,0))-V T (k,x) + Tu(k) . (92) 

Then, defining v{k) := VT(k,<f)^(k,k ,x ,UJ? kak C)), using (J22J)) P3^1. the following holds along 
trajectories 

v(Jfc + 1) - < -T2a(v(A;)) + T[v x + u(k)\ 

for all k > k > 0. From this inequality we can see that the following hold true for any 
k > k Q > 0, 

> a" 1 (x/ x + hf^^l) v(k + l)-v(k)<-Ta(v(k)) (93) 
«(Ar+l)<t;(Ar) +7(^+1^1) . (94) 

Define next ({t) := v(k) + (j^ - k^j (v(k + 1) - v(k)) for all t G [fcT, (Jfc + 1)T), k > k > 0. 

Notice that since ((t) is a linear interpolation from v(k) to u(& + 1), which are always non- 
negative, we have that 

< CO) < max{v(k), v(k+ 1)} for any t G [kT, (k + 1)T ), k > k Q > . (95) 
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Moreover, since trajectories of the system ([18)1. (jl9j) are defined for all k > k > 0, the variable 
C(t) is defined for all t > t > 0. Introduce u(t) := u(k),Vt G [fcT, (A; + 1)T), fc > k Q > 0. Using 
definition of a. T* < 1. (pj jl . (Pf and (pS ) we can write for all t e [fcT, (jfe + 1)T), k > k a > 0: 

C(t) > &{v x + \u [tott) \) v(k) > a' 1 (y x + |wf fco)fc) |) 

=> v(k + 1) - v(k) < —Ta(v(k)) 
C(t) < -Ta(v(k)) 

=► C(*)<-ra(CW), (96) 

where the last inequality follows from the fact that v(k + 1) < w(/c) and hence < v(k) for 
all t e [ifeT, (A; + 1)T), k > k o >0. Using jlSl Corollary 5.1 and 5.2], we have that: 

C(t) < max{/5(Co,t -t ), + |%,*)|)} > 

and using k , x , ^fj^^)) = C(kT), u{kT) = u(k), (3 y and v\, p2|) and the fact that for 

arbitrary a G /C^ we have air + s) < a (2s) + a(2r), Vr, s > we can write 

|0£(fc,fco,a;o,^ o , fe) )| < «r X (max{/3(a 2 (|x |), {k - k )T),afa + |wf fco , fc )|)}) 

< max {a^ 1 (^(a 2 (|a;o|) ) (k - k )T)), a^ 1 o a(2u{) + a^ 1 o a(|a;* fcojfc )|)} 

< max{^(|x |,(fc-A; )r),i/} + 7(|^ oifc) |) , (97) 

for all k > K > 0. 

6 Conclusions 

We have established necessary and sufficient conditions for semiglobal practical asymptotic 
stability of time-varying parameterized cascades. Parameterized systems arise naturally when 
an approximate discrete-time model of a sampled-data plant is used for controller design. Our 
results provide the controller designer with a range of tools that can be used for a systematic 
digital controller design based on approximate discrete-time models. The utility of our results 
was illustrated with a case study, where we have obtained a controller that performs better 
than the emulated (discretized) continuous-time controller designed in jSHj and implemented 
using a sampler and zero-order-hold. 



References 

[1] D. Angeli. Intrinsic robustness of global asymptotic stability. Syst. & Contr. Letters, 
38:297-307, 1999. 

[2] M. Arcak, D. Angeli, and E. Sontag. A unifying integral ISS framework for stability of 
nonlinear cascades. SI AM J. on Contr. and Opt, 40:888-1904, 2002. 

[3] K. J. Astromand Bohn. Numerical identification of linear dynamic systems from normal 
operating records. In P. H. Hammond, editor, Proc. of the 2nd IFAC Symp. on Theory 
of Self-adaptive Control Systems, pages 96-111, Nat. Phys. Lab., Teddington, England, 
1965. 

[4] D. Dochain and G. Bastin. Adaptive identification and control algorithms for nonlinear 
bacterial growth systems. Automatica, 20:621-634, 1984. 



24 



[5] R. Ferretti. High order approximations of linear control systems via Runge-Kutta schemes. 
Computing, 58:351-364, 1997. 

[6] G. C. Goodwin, B. Mclnnis, and R. S. Long. Adaptive control algorithm for waste water 
treatment and pH neutralization. Optimal Contr. Applic. Meth., 3:443-459, 1982. 

[7] L. Griine. Asymptotic behaviour of dynamical and control systems under perturbation and 
discretization. Springer Verlag, Berlin, 2002. 

[8] L. Griine and P. E. Kloeden. Higher order numerical schemes for amnely controlled non- 
linear systems. Numer. Math., 89:669-690, 2001. 

[9] L. Griine and D. Nesic . Optimization based stabilization of sampled-data nonlinear sys- 
tems via their approximate discrete-time models. SIAM, 2002. Accepted for publication 
subject to minor revision. 

[10] A. Isidori. Nonlinear control systems. Springer Verlag, 3rd. edition, 1995. 

[11] Z. P. Jiang. Iterative design of time- varying stabilizers for multi-input systems in chained 
form. Syst. & Contr. Letters, 28:255-262, 1996. 

[12] Z-P. Jiang and Y. Wang. Input-to-state stability for discrete-time nonlinear systems. 
Automatica, 37(6):857-869, 2001. 

[13] Y. Kanayama, Y. Kimura, F. Miyazaki, and T. Naguchi. A stable traking control scheme 
for an autonomous vehicle. In Proc. IEEE Conf. Robotics Automat., pages 384-389, 1990. 

[14] C. Kellett. Advances in converse and control Lyapunov function. PhD thesis, University 
of California, Santa Barbara, CA, USA, June 2002. 

[15] H. Khalil. Nonlinear systems. Macmillan Publishing Co., 2nd ed., New York, 1996. 

[16] I. Kolmanovsky and H. McClamroch. Developments in nonholonomic control problems. 
Control systems magazine, pages 20-36, Dec. 1995. 

[17] D. S. Laila and D. Nesic . Changing supply rates for input-output to state stable discrete- 
time nonlinear systems with applications. Automatica, 2001. Accepted subject to minor 
revision. 

[18] D. S. Laila and D. Nesic . Lyapunov based small-gain theorem for parameterized discrete- 
time interconnected ISS systems. IEEE Trans, on Automat. Contr., 2002. Submitted. 

[19] D. S. Laila, D. Nesic, and A. R. Teel. Open and closed loop dissipation inequalities under 
sampling and controller emulation. European J. of Contr., 8:109-125, 2002. 

[20] A. A. J. Lefeber. Tracking control of nonlinear mechanical systems. PhD thesis, University 
of Twente, Enschede, The Netherlands, 2000. 

[21] W. Lin and C. Byrnes. Passivity and absolute stabilization of a class of discrete-time 
nonlinear systems. Automatica, 31(2):263-267, 1995. 

[22] A. Lorfa. Cascaded nonlinear time-varying systems: analysis and design. Technical report, 
LAG UMR CNRS 5528, 2000. Lecture notes, minicourse at IPN, Mexico, Nov. 2001. 
Available upon request (loria@lss.supelec.fr). 



25 



A. Lorfa and D. Nesic . On uniform boundedness of parameterized discrete-time systems 
with decaying inputs: applications to cascades. Syst. & Contr. Letters, 49(3):163-174, 
2003. In press. 

A. Madani, S. Monaco, and D. Normand-Cyrot. Adaptive stabilization of discrete-time 
strict feedback systems. In Proc. 5th. European Contr. Conf., Karlsruhe, Germany, 1999. 

I. M. Y. Mareels, H. B. Penfold, and R. J. Evans. Controlling nonlinear time-varying 
systems via euler approximations. Automatica, pages 681-696, 1992. 

F. Mazenc and L. Praly. Adding integrators, saturated controls and global asymptotic 
stabilization of feedforward systems. IEEE Trans, on Automat. Contr., 41(11):1559-1579, 
1996. 

S. Monaco and D. Normand-Cyrot. On the sampling of a linear analytic control system. 
In Proc. 24th. IEEE Conf. Decision Contr., pages 1457-1462, Fort Lauderdale, FL, USA, 
1985. 

S. Monaco and D. Normand-Cyrot. On nonlinear digital control, volume 3 of Control, 
chapter in Nonlinear Systems. Chapman & Hall, (A.J. Fossard, D. Normand-Cyrot, Eds.), 
1995. 

P. Morin and C. Samson. A characterization of the Lie algebra rank condition by transverse 
periodic functions. SI AM J. on Contr. and Opt., 40(4): 1227-1249, 2002. 

D. Nesic and D. Angeli. Integral versions of ISS for sampled-data nonlinear systems via 
their approximate discrete-time models. IEEE Trans, on Automat. Contr., 2002. Accepted 
for publication. 

D. Nesic and D. S. Laila. A note on Input-to-State Stabilisation for nonlinear sampled-data 
systems. IEEE Trans, on Automat. Contr., 47(7):1153-1158, 2002. 

D. Nesic and A. Teel. Changing supply functions in input-to-state stable systems: the 
discrete-time case. IEEE Trans, on Automat. Contr., 46(6):960-962, 2001. 

D. Nesic, A. Teel, and P. Kokotovic. Sufficient conditions for stabilization of sampled-data 
nonlinear systems via discrete-time approximations. Syst. & Contr. Letters, 38:259-270, 
1999. 

D. Nesic , A. Teel, and E. Sontag. Formulas relating JCC stability estimates of discrete-time 
and sampled-data nonlinear systems. Syst. & Contr. Letters, 38:49-60, 1999. 

D. Nesic and A. R. Teel. A framework for stabilization of nonlinear sampled-data systems 
based on their approximate discrete-time models. IEEE Trans, on Automat. Contr., 2002. 
Accepted subject to minor modifications. 

E. Panteley, E. Lefeber, A. Lorfa and H. Nijmeijer. Exponential tracking of a mobile 
car using a cascaded approach. In IFAC Workshop on Motion Control, pages 221-226, 
Grenoble, France, 1998. 

E. Panteley and A. Loria . On global uniform asymptotic stability of non linear time- 
varying non autonomous systems in cascade. Syst. & Contr. Letters, 33(2):131-138, 1998. 



26 



E. Panteley and A. Lorfa . Growth rate conditions for stability of cascaded time- varying 
systems. Automatica, 37(3):453-460, 2001. 

G. Sansone and R. Conti. Nonlinear Differential equations. Pergamon Press, Oxford, 
England, 1964. 

R. Sepulchre, M. Jankovic, and P. Kokotovic. Constructive nonlinear control. Springer 
Verlag, 1997. 

C. Simoes, H. Nijmeijer, and J. Tsinias. Nonsmooth stabilizability and feedback lineariza- 
tion of discrete-time sytems. Int. J. Rob. Nonl. Contr., 6:171-188, 1996. 

E. D. Sontag. Remarks on stabilization and Input-to- State stability. In Proc. 28th. IEEE 
Conf. Decision Contr., pages 1376-1378, Tampa, Fl, 1989. 

E. D. Sontag. Comments on integral variants of ISS. Syst. & Contr. Letters, 34:93-100, 
1998. 

E.D. Sontag and Y. Wang. On characterizations of the input-to- state stability property. 
Syst. & Contr. Letters, 24:351-359, 1995. 

A. M. Stuart and A. R. Humphries. Dynamical systems and numerical methods. Cambridge 
University Press, Cambridge, 1996. 



A Some technical proofs 

Proof of Proposition [TJ 

The result is only proved for USC since the proof of UGC follows the same steps. Let ft x G ICC 
come from SP-UAS. Let (A, v) be given and let T{ > be generated using SP-UAS. Without 
loss of generality assume that A > v. Let Ai := (3 X (A, 0) + 1. Let A 2 > be arbitrary. Let 
(Ai, A 2 ) generate K, T 2 * > using conditions of the Lemma. Let T* := minjT^, T 2 *}. Consider 
arbitrary numbers t] G (0, A), e, L > and define 

. fminie, 1} . 
/i:=mm| ,A 

Consider arbitrary \x Q \ < r], \u z \ < fi, T G (0,T*). Note first that without loss of generality 
we may assume that (3 X satisfies fjB3j) . Consequently, SP-UAS implies that whenever \x a \ < A, 
T G (0,T*) then \<f%,(k, k a , x Q , 0)| < max{(3 x (A, 0), v} < (3 X (A,0) < A : for all k > k a > 0. 

With the goal of showing contradiction suppose that there exists k\ G (k a , k + &l,t) such 

that 4>^(k, k a , x , u* ka < A 1 for all k G [k , k\) and (^{k^, k D , x a , uj* ka fe *^) > A^ Then, we 

can write 



\rAk,k ,x ,u z [k ^ k) ) - (j) x T (k,k o ,x o ,0)\ < KT ^ (1 + KT) 

j=0 



U lko,k)\ 



< [(1 + KT) k ~ k ° - 1] \u z \ 

< [e A ' r ^)-l]|^| 

< [e KL -l]\u z \ 

< min{e, 1} (98) 
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for all k G [k Q , k*{}. This contradicts the assumption that 0t(&*, k ,x , 0Jy ko k * 



> Ai since 



0^(A; 1 , k , x Q , k>[k 0) fc*)) 5! 10^(^,^0,2;, 



< ^(A, 0) + 1 = Ai . 

Hence (j5Bj) holds for all > fc + ^l,t, which completes the proof. ■ 
Proof of Claim 1: 

Let Tj* > come from inequalities and (|48jl. Let /i, L,T 2 * come from the PE condition. 
Let u>m come from the conditions of the proposition. Let T 3 * > and T 4 * > be such that 



l-e- 



<2, VTe(0,T ; 



37) 



1 - e~ T 1 

~^^>2' VTG(0,T 4 *). 



(99) 



Let T 5 * := A ^ e _ e - L) and, finally, define T* := min{T*, T*, T 3 *, T 4 *, T*}. Consider arbitrary x eR 2 , 
k > and T G (0, T*). The first inequality in (fl9*J) holds using our choice of T 3 *: 



i=o 



i=k 



1-e 



-T 



< 2w 



M = : C 3 



(100) 



The second inequality in (jl9*j) follows from the PE condition. In particular, denote kj := k+j£L,T 
and then we can write: 



i=k 



,{k-i) Tu 2 



00 

E r E 

3=0 i=k+ji L , T 



kj+e L ,T 



3=0 



> 



> 



E 

j=o 



jL-£ L , T T 



e J e 



T E 



1-e- 



-/i =: c 4 . 



(101) 
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We show next that (|5U|) holds. Using similar calculations as in ( [100)1 . (jl01|) and our choice of 
T* and T* = f we have: 

w T(k + i,ye-^ k x e )-w T (k,y e ) = _ ^ e ^yr <{ye _ TUrkXe? 

i=k+l i=k 

oo 

= <VI - E ^ l)T < l(Ve ~ Tu rk x e f - yj] 

i=k+l 

oo 

-(1 - e- T ) e^-Wuttf 



i=k+l 

oo 



< ^l k yl + ( ^M T E e ^ ) [ 2UjM W + Tw M X e] 

3=0 



1 - 

- 1 -i)T 2 2 



T 

i=fc+l 



^ 2 2 i M fn I II I i rrr 2 2\ ^12 

- W r u Ve + I _ e _ T (2w M ke| |j/ e | + Tw M X e ) - —y e . 

Finally, using (j99|) and completing the squares we obtain that there exist K 2 , a y > such that 

< u 2 rk yl + K 2 x 2 e - a y y 2 e , (102) 



W T (k + l,y e - Ur k x e ) - W T {k, y eJ < ^ 2 2 



which completes the proof of the claim. 
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A.l Proof of Proposition [2] 

First, note that using (J32*|) . (J33)l we can write 

V T (k+l,f T (k,x,z))-V T (k,x) = V T (k + l,f T {k,x,0))-V T (k,x) 

+V T (k + 1, f T (k, x, z)) - V T (k + 1, / T (A;, x, 0)) 
< TU\z\MV T (k,x)) + T%(\z\). 

Assume for the time-being that from this inequality and the conditions of the Proposition, there 
exist at, a 2 G /Coo, c > and for each T G (0, T*) there exists W T ■ K>o x ^ Ux — > ^>o such that 
all x G M n , k > and T G (0, T*) we have that 

«i(|x|) < W" T (M) < a 2 (\x\) + c (103) 
W T (fc + l,/ T (M,z))-W- T (fc,x) < 7>(|z|). (104) 
Then, defining w(k) := Wx(k, (f)^(k, k , x , fc ^)) and using ()104j) we can write 

< io(fc) + T^(|o;f feo)fc) |), Wk>k o >0 
and consequently, by induction it follows that 

k-l 

w(k) < w(k Q ) MK ,fc)|), VA; > A; > . 

Next, using ()103|) we obtain that 

k-l 

ai (|0r(fc,fco,x o ,a;f fcoiJfc) )|) < a 2 (|x |) + c + T ^ //(|wf fcojfc) |), VA: > k a > , (105) 

t = ko 

hence, UGB follows from fl35jl and the fact that a±, a 2 G /Co- 

It is only left to prove that under the conditions of the proposition ()103j) and ()104)) hold. 
To that end let &i, a 2 , (p, 7i, 72 come from Assumption^ Define 

Let p G /Cqo be defined as 

p(s) := / q(r)dT . 



Note that p G /Cqo and p is differentiable with -f(s) = q(s). Also, q(-) is non-increasing 
since ip G /Cqo. Introduce Wt(A:, x) := p(Vr(A;, x)). Then, it is clear that ()103|) holds with 
ai(s) : = p o fii(s), «2(s) = p o 2a 2 (s) and c = p(2c). We now show that (|l()4jl holds. 

We consider two cases. If Vr(k + 1, fr(k, x, z)) < Vr(k, x), then since p G /Coo, we have: 

W T (k + 1, x, z)) - W T (k, x)<0<T(j l (\z\) + j 2 (\z\)-^) • (106) 



On the other hand, if Vr(fc + 1, /r(&, x, z)) > Vr(k, x), then using the Mean Value Theorem we 
can write 

W T (k+l,f T (k,x,z))-W T (k,x) = q(V*)[V T (k+l,f T (k,x,z))-V T (k,x)] 

< Tq(V*)M\z\MV T (k,x))+%(\z\)} (107) 



30 



where Vr(k, x) < V* < Vx(k + 1, fr(k, x, z)). Since q is a non increasing function, we have 

q(V T (k,x)))>q(V*) , 

and hence 

W T (k + l,f T (k,x,z))-W T (k,x) < Tq(V T (k,x))[U\z\MV T (k,x))+U\z\)} ■ 
Now, if Vr{k,x) < 1 then we have that 

M\*\Mv T (k,x)) + j2(\z\) 



W T (k+l,f T (k,x,z))-W T (k,x) < T 

< T 



7i(klMl)+72(kl) 
¥>(1) 



= T (ti(N) +72(1^1)^) • (108) 

Otherwise, if Vt(&, x) > 1 then, 

,7i(klMMM)) + 72(1*1) 



^ T (fc + i,/ T (fc,x,^))-iy r (fc,x) < t- 



<p{Vr{k,x)) 
U\A) 



<p(V T (k,x)) 



r(7i(N) + 7 2 (N)^j] • (109) 



The proof is completed by combining (jlOfjjl . (|l()8jl and (|l()9jl . 
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